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Abstract 


Tliip  project  is  based  on  a  program  that  uses  FORTRAN  77  to  aii- 
rlyee  resistance  versus  temperature  data.  Tlie  data  aie  fit  by  two  lines 
vitli  dilTereiit  slopes  and  intercepts.  A  toial  for  tlie  fits  is  obtained. 
^Vllen  tlie  total  is  minimised,  the  critical  lem))erature  is  assigned  to 
the  point  where  the  data  are  broken  into  the  intersecting  two  lines.  The 
program  also  informs  the  user  of  the  temperature  at  which  those  lines 
i  itersect.  A  sensitivity  study  was  performed  to  find  the  relationship 
I  etweeii  uncertainties  in  the  resistance  and  critical  temperature.  This 
c-^nalysis  showed  the  resistance  needed  to  be  measured  with  less  than 
1%  error  to  allow  the  T,-  to  be  found  with  an  error  of  three  Kelvins. 


1 


Introduction 


The  recent  discovery  of  high  temperature  superconductors  has  led  to  the 
study  of  how  fhese  superconductors  work  [1],  Of  particular  interest  is  the 
critical  tempei  aturc  (Tc)  and  slope  of  the  resistance  versus  temperature 
graph  below  tie  critical  temperature.  This  paper  reports  on  a  computer 
program  that  .its  two  lines  to  resistance  versus  temperature  data  from  any 
superconductor.  The  program  identifies  the  critical  temperature  and  slope 
below  the  critical  temperature  of  the  superconductor. 

Another  program  was  made  to  provide  sample  data  sets  that  have  noise 
added  to  them.  The  noise  was  provided  by  a  random  number  generator 
that  sampled  8  gaussian  distribution  with  zero  mean  and  unit  variance.  The 
program  output  was  used  to  study  the  relationship  between  the  uncertainties 
in  the  resistan  :e  values  and  uncertainties  in  the  critical  temperature. 

2  DESCRIPTION  OF  PROGRAM 

This  section  describes  the  aspects  of  the  programs,  and  details  the  theory 
behind  the  code.  The  listings  have  been  labeled  to  facilitate  following  the 
description.  T  le  program  that  fits  the  lines  (THIN. FOR)  is  addressed  first. 

In  the  main  program,  part  1,  the  variables  are  defined  and  the  data  and 
output  files  are  opened.  The  files  are  briefly  described  in  the  comment  block 
below  the  statements.  The  variables  C  and  CH BCK  are  used  in  the  least- 
squares  subroutine  and  are  addressed  later.  Part  2  counts  the  number  of 
data  points  in  the  data  set.  It  also  asks  the  user  if  there  is  an  error  value 
for  each  data  point.  Part  3  calls  the  subroutine  LSTSQR  and  checks  to 
determine  if  another  data  break  is  needed. 

Part  4  of  t’le  program  is  the  subroutine  LSTSQR  (least-squares  fit)  [2]. 
This  subroutine  fits  two  lines  with  different  slopes  and  intercepts  to  the  data. 
It  also  determines  how  closely  the  lines  fit  the  data  by  finding  the  total  x* 
of  the  two  lines.  The  x*  is  the  residual  value  divided  by  the  total  standard 
deviation  for  eJl  the  points.  The  smaller  x*,  the  closer  the  fitted  line  is  to 
the  actual  data.  This  least-  squares  fit  is  repeated  by  continually  moving 
the  point  at  wiiich  the  data  are  broken  into  the  two  lines  back  towards  zero 
temperature.  The  Tc  is  assigned  by  the  user  when  the  x^  is  minimized  for 
both  lines. 

In  Part  4a,  two  counters,  LINENUM  and  COUNT,  are  initialized  to 
one.  LINENUM  simply  denotes  which  line  is  being  fit,  and  COUNT  is  the 
current  number  of  the  data  point  being  read  in.  The  R  and  T  data  are  read 


in  from  file  RVST2.DAT.  The  temperature  data  is  measured  in  Kelvins  and 
the  resistance  data  is  a  relative  resistance  with  zero  denoting  no  resistance 
and  ore  denoting  the  maximum  resistance  for  a  given  material.  The  data 
are  an  anged  in  two  columns  with  the  independent  variable  (temperature) 
on  the  left  and  the  dependent  variable  (resistance)  value  on  the  right.  The 
columns  are  separated  by  a  comma.  The  values  in  RVST2.DAT  can  be 
integei  or  real.  The  data  are  read  into  two  arrays  A'  and  T.  The  values  of 
A  and  Y  are  summed,  and  this  information  is  used  in  determining  the  slope 
and  intercept  for  the  line  that  is  being  fit. 

Pa-t  4b  removes  the  first  data  point  in  each  set  until  a  second  point  is 
input.  This  fills  the  requirement,  that  a  line  needs  two  points  to  be  defined. 
Part  4  :  determines  the  standard  deviation  of  the  data  point.  The  standard 
deviation  can  be  user  defined,  or  can  be  approximated  from  the  data  by 
taking  the  square  root  of  the  data  point.  The  standard  deviation  is  later 
used  to  determine  the  ■  Part  4d  checks  to  see  if  there  are  more  data  points 
to  be  f  t.  If  not,  the  subroutine  skips  to  the  end.  The  slope  and  intercept  are 
deterrr  ined  in  4e,  and  the  values  are  written  to  the  external  file  OUT.OUT 
in  the  format  prescribed  in  line  55.  Next,  the  residual  value  is  computed  in 
part  4f. 

Pa  ’t  4g  is  the  test  to  determine  if  one  line  has  been  fit.  Inside  the  sub¬ 
routine,  the  program  checks  to  see  if  a  dummy  variable,  C  (a  counter),  is 
greater  than  the  predetermined  number  of  checks  desired.  The  number  de¬ 
sired  in  set  in  the  main  program  (CHECK).  If  there  hats  not  been  enough 
data  p  Dints  run  through  the  subroutine,  then  the  program  reads  in  another. 
If  that  is  the  point  at  which  the  user  wants  the  first  line  break,  then  the  sub- 
routin«!  assigns  the  temperature  value  to  that  of  a  variable  called  the  critical 
tempei  ature  and  begins  to  fit  the  second  line.  The  same  process  is  followed 
for  the  second  line,  and  a  new  slope,  intercept,  and  are  determined. 

A  lotal  X*  is  defined  as  the  root  mean  sum  of  the  individual  x*  of  the 
lines.  This  value  is  written  to  OUT.OUT,  and  the  process  is  repeated  with 
the  br'iak  in  the  lines  being  moved  back  one  data  point  from  the  previous 
run.  The  user  seeks  to  minimize  the  total  x^i  and  thus  pick  the  proper 
place  for  the  T^.  The  program  outputs  the  slope  and  intercept  for  both 
lines.  These  actions  are  incorporated  in  parts  4g,  4i,  and  4j. 

Pa:  t  4h  is  inserted  to  determine  the  errors  associated  with  the  resistance, 
the  critical  temperature,  and  the  slope  of  the  lines.  All  these  values  are 
output  to  different  files  as  prescribed  at  the  top  of  the  main  program.  The 
error  i:i  the  resistance  is  read  from  an  external  file  (  S1GMAR.DAT)  as  a 
percen  tage  and  then  multiplied  by  the  the  value  of  R.  Pait  4K  merely  puts 


the  output  data  in  a  form  that  is  aesthetically  pleasing. 


2.1  NOISID 

The  second  program,  NOISE,  generates  random  numbers  based  on  the  gaus- 
sian  distribution  of  the  data  points  and  adds  random  ’’noise”  to  the  data 
points  [3|.  These  points  can  then  be  read  by  the  first  program  as  input. 
Tins  allows  THIN. FOR  to  be  tested  on  data  that  is  corrupted  by  noise  of 
varying  amour  ts. 

In  part  1  tne  data  are  read  into  the  program  exactly  as  in  THIN. FOR. 
Part  2  either  adds  or  subtracts  "noise”  from  the  resistance  values.  The 
choice  of  addition  or  subtraction  is  random.  Part  3  randomly  changes  the 
seed  from  which  the  random  number  is  generated. 

There  are  two  subroutines,  URAND  (part  5)  and  GAUSSN  (part  4). 
URANU  generates  a  random  number  between  zero  and  one.  It  does  this  by 
a  series  of  divisions  and  multiplications.  The  remainders  of  the  iivisions  are 
used  to  produce  the  seed  which  is  used  with  the  intrinsic  command  FLOAT 
to  produce  a  random  number.  This  number  is  in  turn  used  by  GAUSSN 
to  produce  a  normally  distributed  random  number  with  unit  variance.  The 
program  listin.^s  are  found  in  Appendices  A  and  B. 


3  RESULTS 

The  sample  dr.ta  in  RVST2.DAT,  found  in  Appendix  D,  are  characteristic 
of  R  versus  T  graphs  [4].  The  data  have  been  hand  digitized  from  the  orig¬ 
inal  paper  which  sustituted  preiseodyium  into  a  bismuth  strontium  calcium 
copper  oxide  superconductor.  The  resistance  versus  temperature  output  in 
Appendix  C  dictates  that  the  line  be  broken  at  approximately  104  deg  K. 
This  is  very  close  to  the  middle  of  the  "elbow”  in  the  actual  data  between 
llSand  97K  as  seen  graphically  in  Figure  1. 


3.1  UNCERTAINTIES  IN  PARAMATERS 

The  uncertainly  in  the  slope  was  determined  as  a  function  of  the  uncertain¬ 
ties  in  resistance.  Following  Bevington, 


RELATIVE  RESISTANCE 


R  vs.  T: 

Data  and  Fitted  Lines 


Figure  1 


=  -  2iVi,  IZ  +  ( 51  2:.)* 

j=i  ^ 

where  A  =  iV  X)  -  (  Z)  *i)* 


with  AT  =  16,  ffr  =  0.003,  and  A  =  0.64,  we  obtain  jS]. 

Th2  uncertainty  in  the  critical  tempertitnre  (Tr)  was  found  by  two  meth¬ 
ods.  The  first  was  found  by  using  the  slope  of  the  first  line  and  the  (Tr 
the  po  nt  where  the  lines  were  broken  (point  slope  method)  The  second  was 
found  by  using  the  iut ’rsection  of  the  two  lines  and  incorporating  in  the 
calcule  tion  the  errors  in  the  slopes  of  both  lines  and  their  corresponding  y- 
intercepts  (intercept  method).  The  equations  for  these  methods  are  shown 
below. 


Point  Slope 


M  s  (  he  respective  slope  for  each  line  and  b  is  the  respective  y  intercept. 
Appendix  C  shows  m  and  b  for  each  fitted  line. 


TABLE  1 


SIGMA  R  /  R 

SIGMA  R  @  TC 

SIGMA  TC 
(PT.  SLOPE) 

SIGMA  TC 
(INTERSEPT) 

SIGMA  SLOPE 

0.9 

0.58 

35.48  K 

73.09 

K 

0.014214/K 

0.7 

0.45 

27.60  K 

56.86 

K 

0.011055/K 

0.4 

0.26 

15.77  K 

32.49 

K 

0.006317/K 

0.2 

0.13 

7.88  K 

16.23 

K 

0.003159/K 

0.1 

0.06 

3.94  K 

8.12 

K 

0.001580/K 

0.05 

0.03 

1.97  K 

4.06 

K 

0.000790/K 

SIGMA  SLOPE  VS.  SIGIIA  R 
m  =  0.02256/K 

SIGMA  TC  VS.  SIGMA  R 
m  =  58.43  K 


SLOPE  OF 
LEFT  LINE 
0.02/K 


TC 

104K 


Sigma  Tc  (K) 


Sigma  Tc  vs.  Sigma  R  /  R 


Point  Slope  Method 

— I — 

Intercept  Method 


Sigma  R  /  R 


Figure  2 


Table  J  shows  <7r..  for  different  values  of  Figure  2  shows  the  same 
data  graphical  y.  Figure  2  shows  a  ^  of  about  2  -  5  of  approximately  3  K. 

4  CONCLUSIONS 

Using  the  hand  digitized  data  in  RVST2.DAT  (4),  the  program  THIN. FOR 
fits  two  lines  to  the  resistance  versus  teinf>erature  data  produced  by  thin 
film  supcrconcuctors.  The  program  determines  the  slope  and  intecepts  of 
the  two  lines  a  id  the  critical  temperature.  The  program  also  allows  the  user 
to  determine  where  one  line  ends  and  the  other  starts,  as  appropriate  with 
the  user’s  neecs. 

The  uncerlanty  in  the  critical  temperature  was  found  to  be  related  to 
the  relative  uncertainty  in  the  resistance.  It  was  found  that  a  few  percent 
relative  measu  "ement  of  resistance  is  necessary  to  determine  Tc  with  an  error 
of  3  Kelvins.  It  was  also  found  that  the  percent  error  in  the  slope  is  small 
compared  to  u  r  and  . 
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APPENDIX  A 


**********************************************  ^*  *******  *  *************** 

* 

THIN  FILMS  PROGRAM 
BRIAN  A.  DOYLE,  CS-10,  X4500 
29  NOV  89 

This  program  analyzes  resistance  versus  temperature. 

It  can  be  used  in  conjunction  with  superconductors. 

A  step  by  step  analysis  of  the  program  can  be  found 
in  the  "Description  Of  Program"  section  of  the  paper. 

* 

★ 

****★**★★★*★***★★★★★*★***** *★★★★**********★*★**★*★★*****★*★★★★**★**★***★ 
PROGRAM  THIN 

IMPLICIT  LOGICAL  (A-Z) 

INTEGER  I,N,  COUNT,  NUMBERDAT,  DEV 

REAL*8  X(500),  Y(500),  SLOPE,  YINT,  YNEW,  RESID,  SUMX,  SUMY 

REAL*8  SUMXY,  SUMXX,  SUMRESID,  CHITOT,  CHECK,  C 

DATA  COUNT,  SUMX,  SUMY,  SUMXY,  SUMXX,  SUMRESID  /I, 5*0.0/ 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


TITLE: 
AUTHOR : 
DATE: 
OVERVIEW: 


* 

* 

* 

* 

* 

* 

* 

* 


OPEN  (UNIT=10,FILE=  'RVST.DAT',  STATUS  =  'OLD') 
OPEN  (UNIT=11,FILE=  'OUT. OUT',  STATUS  =  'NEW') 

OPEN  (UNIT=12,FILE*  'FIT. DAT',  STATUS  =  'NEW') 

OPEN  (UNIT=13,FILE=  'ERR.DAT' ,  STATUS  *  'OLD') 

OPEN  (UNIT=14,FILE*  'RVST 2 .DAT' .STATUS  «  'OLD') 
OPEN  (UNIT=15,FILE=  'ERRSLOP. OUT' .STATUS  =  'NEW') 
OPEN  (UNIT=16,FILE=  ' ERRTC . OUT ', STATUS  =  'NEW') 
OPEN  (UNIT=17,FILE=  'SIGMAR.DAT' .STATUS  =  'OLD') 

10  =  data  points  used  to  count 

11  =  output  of  thin 

12  =  data  points  of  the  fitted  lines 

13  =  error  for  each  data  point 

14  =  data  points  used  for  least  squares  fit 

15  =  error  associated  with  the  slope 

16  =  error  associated  with  the  Tc 

17  =  percent  error  associated  with  R  (it  is  altered) 
C  =  0.1 

CHECK  =1.32 


1 


DO  9  I  =  1,1000 

READ  (10,*)  X(I) 

IF  (X(I)  .GE.  0)  THEN 

NUMBERDAT  =  NUMBERDAT  +  1 

ELSE  2 

write(*,*)  'There  are ' , nunberdat , '  data  points.' 

CLOSE  (10) 

GOTO  10 
ENDIF 

9  CONTINUE 


Switch  to  enter  error  or  use  square  root  of  data  point  * 

WRITE (*,*)  'Input  (1)  if  you  have  an  error  value  for  each  data  poi 
&nt . ' 

WRITE(*,*)  'Otherwise,  enter  (2)' 


*  00 


READ  (*,*)  DEV 

*  Go  thru  process  of  fitting  lines  until  there  are  no  more  data  points 
1  =  0 

DO  8  I  =  1,  NUMBERDAT 

1  CALL  LSTSQR(X,  Y,  NUMBERDAT,  CHITOT,  CHECK,  C,  COUNT,  DEV) 

IF  (CHITOT  .GT.  C)  THEN 

CHECK  =  CHECK  -.091  3 

C  =  CHITOT 
COUNT  =  1 
ENDIF 
CONTINUE 

Closing  the  files  previously  opened. 

CLOSE  ( 17 ) 

CLOSE  (16) 

CLOSE  (15) 

CLOSE  (14) 

CLOSE  (13) 

CLOSE  (12) 

CLOSE  (11) 

CLOSE  (10) 

STOP  'OUTSTANDING' 

END 


LSTSQR 

BRIAN  A.  DOYLE,  CS-10 ,  X4500 
29  NOV  89 

This  subroutine  uses  the  least  squares  fit  method  to  fit 
the  R  and  T  data,  in  an  external  file,  to  a  line. 

See  pages  382-388,  453  in  Structured  Fortran  by  Etter. 

SUBROUTINE  LSTSQR ( X , Y , NUMBERDAT , CHITOT , CHECK , C , COUNT , DEV ) 


it 

* 

TITLE 

* 

AUTHOR 

* 

DATE 

it 

OVERVIEW 

* 

* 

SOURCES 

IMPLICIT  LOGICAL  (A-Z) 

INTEGER  COUNT,  N,  I,  LINENUM,  NUMDATA,  NUMBERDAT,  DEV,  RUFF,  ARF 
INTEGER  GOODFIT,  Q 

REAL*8  X(500),  Y(500),  SLOPE,  YINT,  YNEW,  RESID,  SUMX,  SUMY 
REAL*8  SUMXY,  SUMXX,  SUMRESID,  CHIl ,  SIGMA(500),  SIGMATOT,  CHITOT 
REAL*8  CHI2,  CHECK,  C,  SIGTC,  SIGSLOPE,  ERRY,  SIG,  DELTA,  SUMX2 
REAL*8  SLOl,  SL02,  YINTl ,  YINT2 ,  INTSEC 


GOODFIT  =  0 
LINENUM  =  1 
COUNT  =  1 

Reading  in  data  points.  ** 

This  is  so  I  can  mess  with  NUMDATA  and  still  know  how  many 
points  there  are. 

NUMDATA  =  NUMBERDAT 


5  READ  (14,*)  X(COLT<T),  Y  (COUNT) 

SUMX  =  SUMX  +  X( COUNT) 

SUMY  =  SLUY  Y(  COUNT)  4a 

SUMXY  =  SUMXY  +  X( COUNT)  *  Y( COUNT) 

SUMXX  =  SUMXX  +  X( COUNT)  *  X( COUNT) 

*♦  This  is  to  t.hrow  out  the  first  data  point  of  each  line  because 
•  it  produces  bad  results. 

IF  (COUNT  .EQ.  1)  THEN 


COUNT  =  COUNT  +  1 

GOTO  5  4b 

END  IF 

COUNT  =  COUNT  +  1 
N  =  COUNT  -  1 

Figuring  the  standard  deviation  for  each  point.  ** 

The  SD  is  either  read  in  or  the  square  root  is  used. 

IF  (DEV  .EQ.  1)  THEN 
READ  (13,*)  SIGMA(N) 

ELSE 

SIGMA(N)  =  DSQRT(X(N)) 

SIGMATOT  =  SIGMA(N)  +  SIGMATOT 
ENDIF 

If  there  are  no  more  data  points,  go  to  the  end. 

IF  (COUNT  .GT.  NUMDATA)  THEN 
GOTO  53 
ENDIF 

IF  (N  .LE.  1)  THEN 
N  =  2 
ENDIF 

Determining  slope  and  y  intercept  of  one  line.  ** 

SLOPE  *  (SUMX  ★  SUMY  -  REAL(N)  *  SUMXY)  / 

&  (SUMX  *  SUMX  -  REAL(N)  *  SUMXX) 

YINT  *  (SUMY  -  SLOPE  *  SUMX)  /  REAL(N) 

IF  (LINENUM  .EQ.  1)  THEN 
SLOl  =  SLOPE 

YINTl  *  YINT  4e 

ELSE 

SL02  =  SLOPE 
YINT2  =  YINT 
ENDIF 

WRITE (11,*)  'THE  LINEAR  EQUATION  IS  FOR  LINE  #', LINENUM 
WRITE (11, 55)  SLOPE,  YINT 

Calculating  a  predicted  value  for  y  and  determining  how  good  the 
fit  is,  based  on  this  value. 

DO  65  I  =  1,N 

YNEW  =  SLOPE  *  X(I)  +  YINT 
IF  (C  .GE.  CHECK)  THEN 
WRITE  (12,91)  X(I),  YNEW 
ENDIF 

RUFF  =*  ARF  +  I  -  l.ODO 

IF  (LINENUM  .EQ.  2  .AND.  COUNT  .EQ.  NUMDATA)  THEN 

WRITE  (12,92)  X(RUFF),  YNEW  4f 

ENDIF 

RESID  =  Y(I)  -  YNEW 

SUMRESID  =  SUMRESID  +  RESID  *  RESID 
WRITE(11,60)  X(I),  Y(I),  YNEW,  RESID 
CONTINUE 

WRITE  (11,*) 

WRITE  (11,70)  SUMRESID 

Checking  to  see  if  the  fit  is  good  for  one  line  or  if  it  is 
trying  to  use  a  point  on  the  second  line.  If  it  is  not  good, 
the  fit  for  a  second  line  is  started. 

IF  (C  .LT.  CHECK)  THEN 


4C 


4d 


C  =  C  +  .1 
GOTO  5 

ELSE 

When  one  line  has  been  fit,  the  following  hapoens. 

C  =  -5.0 

WRITE(*,*)  'AT  DATA  POINT  NUMBER  ',N,',  ONE  LINE  HAS  BEEN  FIT 
WRITE(11,75)  X(N) 

CHIl  =  DSQRT( (SUMRESID  /  SIGMATOT)  **  2)  4g 

CHIl  =  CHIl  *  100. ODO 
WRITE  (11,80)  CHIl 
CHIl  =  CHIl  /  100. ODO 

For  t.he  best  fit  (goodfit),  the  errors  are  written  to  another  fi] 
GOODFIT  =10 

IF  (N  .EQ.  GOODFIT)  THEN 
READ  (17,*)  ERRY 
Q  =  N  -  2 

write(*,*)  'y(Q)  =  ',y(Q) 

SIG  =  ERRY  *  Y(Q)  4h 

SIGTC  =  SIG  /  SLOPE 
This  coroes  from  page  114  in  Bevington 
SUMX2  =  SUMX  **  2 
DELTA  =  (N  *  SUMXX)  -  SUMX2 
SIGSLOPE  =  DSQRT(N  *  (SIG  **  2)  /  DELTA) 

WRITE (15, 93)  SIG,  SIGSLOPE 
WRITE (16, 91)  SIG,  SIGTC 
END  IF 

NUMDATA  =  NUMDATA  -  N 
SUMX  =  O.ODO 
SUMY  =  O.ODO 

SUMXY  »  O.ODO  4i 

SUMXX  =  O.ODO 
SUMRESID  =  O.ODO 
LINENUM  =  2 
ARF  =  COUNT 
COUNT  =  1 
GOTO  5 
ENDIF 

Calculating  chi  for  the  2nd  line  and  the  total  chi  squared. 

IF  (LINENUM  .EQ.  2)  THEN 

CHI2  =  DSQRT( (SUMRESID  /  SIGMATOT)  **  2) 

CHITOT  =  DSQRT(CHI1  **  2  +  CHI2  **  2) 

CHI2  =  CHI2  *  100. ODO 
WRITE  (11,85)  CHI2 
CHI2  =  CHI2  /  100. ODO 
WRITE  (11,87) 

CHITOT  =  CHITOT  *  100. ODO 
WRITE  (11,88)  CHITOT 
CHITOT  =  CHITOT  /  100. ODO 

WRITE  (11,89)  4j 

ENDIF 

SUMX  =  O.ODO 
SUMY  =  O.ODO 
SUMXY  =  O.ODO 
SUMXX  *  O.ODO 
SUMRESID  =  O.ODO 
WRITE  (12,*)  '*★★***★★★' 

INTSEC  »  (YINT2  -  YINTl )  /  (SLOl  -  SL02) 

WRITE  (*,*)  'THE  INTERSECTION  IS  AT  ', INTSEC,'  K' 


54 

^^5 


RETURN 

This  is  only  to  make 
FORMAT  (IX, 'Y  = 


'  ,F6 


it  look 
2  ,  '  X 


nice. 


,F6.2) 


60 

FORMAT 

(IX, 

F6.2,6X,F6.2,6X,F6.2,7X,F6.2) 

■;  70 

FORMAT 

(IX, 

'RESIDUAL  SUM  =' ,  F6 . 2 ) 

■  75 

FORMAT 

(IX, 

'THE  CRITICAL  TEMPERATURE  IS  ', 

F6. 

2  , '  DEGREES 

80 

FORMAT 

(IX, 

'THE  VALUE  OF  CHI  FOR  1ST  LINE: 

t 

f 

F7.4,  /) 

*:  85 

FORMAT 

(IX, 

'THE  VALUE  OF  CHI  FOR  2ND  LINE: 

/ 

/ 

F7.4) 

87 

FORMAT 

(IX, 

/***■*★*★★**★*****★****★★★★**** 

') 

88 

FORMAT 

(IX, 

'TOTAL  CHI  SQUARED  FOR  BOTH  LI 

NES: 

',  F7.4) 

89 

FORMAT 

(IX, 

^ ★★*★★***★*★*»*★**★*****★★**** 

91 

FORMAT 

(5X, 

F6.2,  3X,  F6.2) 

92 

FORMAT 

(5X, 

F6.2,  3X,  F6.2) 

93 

FORMAT 

(5X, 

F6.2,  3X,  F8.6) 

K') 


4k 


APPENDIX  B 


it 

★ 

*  TITLE:  RANDOM  NOISE  GENERATOR 

★ 

*  AUTHOR:  DEPARTMENT  OF  ASTRONAUTICS,  USAFA 

★ 

*  DATE:  16  JAN  90 

it 

*  OVERVIEW:  This  program  generates  random  numbers  by  sampling 

*  a  gaussian  distribution  with  zero  mean  and  unit 

*  variance.  The  program  adds  or  substracts  values 

*  from  R  vs.  T  data  to  create  "noise”  for  a  sensitivity 

*  check. 

★★♦★★★★★★★★♦★★★★★★★★★★★★★★★★★★♦★★★★★It*********************************** 


PROGRAM  NOISE 

REAL*8  YNOISE(500) ,X(500) ,Y(500) ,TEMP,TEMP2 
INTEGER  NUMBERDAT , SEED , I , COUNT ,  ADDOR 

OPEN  (UNIT  =  5, FILE='RVST. DATS  STATUS  =  'OLD") 

OPEN  (UNIT  =  6,  FILE  =  'RAND. DAT",  STATUS  =  'NEW') 

ADDOR  =  1 
NUMBERDAT  =  1 
DO  9  I  *  1,1000 

READ  (5,*)  X(I) 

IF  (X(I)  .GE.  0)  THEN 

NUMBERDAT  =  NUMBERDAT  +  1 

ELSE  1 

NUMBERDAT  =  NUMBERDAT  -  1 
CLOSE  ( 5 ) 

GOTO  10 
END  IF 

9  CONTINUE 

10  OPEN  (UNIT=5,  FILE  =  'RVST.DAT',  STATUS  =  'OLD') 

SEED  ==  824064364 
DO  11  COUNT  =  1, NUMBERDAT 

READ  ( 5 , * )  X ( COUNT ) ,  Y ( COUNT ) 

SIG  =  Y(COUNT) 

TEMP  =  URAND(SEED) 

TEMP2  =  GAUSSN(SIG,  SEED)  /  100 

IF  (ADDOR  .EQ.  2)  THEN 

Y( COUNT)  «  Y( COUNT)  +  TEMP 2  2 

ELSE 

Y( COUNT)  =  Y( COUNT)  -  TEMP2 
ENDIF 

WRITE  (6,50)  X(COUNT),  Y(COUNT) 

IF  (TEMP  .GT.  .5)  THEN 
ADDOR  =  2 

SEED  =  SEED  +  (1000  *  TEMP) 

ELSE  3 

SEED  =  SEED  -  (1000  ♦  TEMP) 

ENDIF 

11  CONTINUE 


CL0SE(5) 

CL0SE(6) 

FORMAT(10X,  F7.2,  F7 . 2 ) 

STOP  'DONE' 

END 

REAL  FUNCTION  GAUSSN(SIG,  SEED) 

IMPLICIT  LOGICAL  (A-2) 

REAL*8  GNOIZ,  SIG 
INTEGER* 4  SEED,  I 

SIG  =5.0 
GNOIZ  =  0. 

DO  10  I  =  1,12  4 

GNOIZ  =  GNOIZ  +  GRAND (SEED) 

CONTINUE 

GAUSSN  =  SIG  *  (GNOIZ  -  6.0) 

RETURN 

END 

REAL  FUNCTION  GRAND (SEED) 

IMPLICIT  LOGICAL  (A-Z) 

INTEGER* 4  B2E15 , B2E16 ,MODLUS ,HIGH15 ,HIGH31 ,LOW15 , LOWPRD 
INTEGER*  4  MULT 1 , MULT 2 , OVFLOW , SEED 
DATA  MULTI, MULT2/24112, 26143/ 

DATA  B2E15,B2E16,MODLUS/32768, 65536, 2147483647/ 

HIGH15  =  SEED/B2E16 

LOWPRD  =  (SEED  -  HIGH15*B2E16 )  *  MULTI 
LOW15  =  LOWPRD/B2E16 
HIGH31  =  HIGH15  *  MULTI  +  LOW15 
OVFLOW  =  HIGH31/B2E15 

SEED  =  (((LOWPRD  -  LOW15  *  B2E16)  -  MODLUS)  +  5 

&  (HIGH31  -  OVFLOW  *  B2E15)  *  B2E16)  +  OVFLOW 

IF  (SEED  .LT.  0)  SEED  =  SEED  +  MODLUS 

HIGH15  *  SEED/B2E16 

LOWPRD  =  (SEED  -  HIGH15  *  B2E16)  *  MULT2 
LOW15  =  LOWPRD  /  B2E16 
HIGH31  =  HIGH15  *  MULT2  +  LOW15 
OVFLOW  =  HIGH31  /  B2E15 

SEED  =  (((LOWPRD  -  LOW15  *  B2E16)  -  MODLUS)  + 

&  (HIGH31  -  OVFLOW  *  B2E15)  *  B2E16)  +  OVFLOW 

IF  (SEED.LT.O)  SEED  =  SEED  +  MODLUS 

GRAND  =  FLOAT ( 2* (SEED/256)  +  1 ) /16777216 . 0 

RETURN 

END 


APPENDIX  C 


This  appendix  contains  the  output  of  the  program 
THIN. FOR  for  the  resistance  versus  temperature  data 
contained  in  Appendix  D.  The  critical  value  is  the  total 
chi  squared  for  both  lines  which  is  denoted  by  two  rows  of 
stars.  This  value  informs  the  user  of  the  best  least- 
squares  fit  match  and  its  associated  critical  temperature. 
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THE  LINEAR  EQUATION  IS  FOR  LINE  # 

V  =  .03  X  ♦  -I .96 

THE  CRITICAL  TEMPERATURE  IS  90.00  DEGREES 
THE  VALUE  OF  CHI  FOR  1ST  LINE:  .0000 


TOTAL  CHI  SQUARED  FOR  BOTH  LlNESi  .0124 


N 


UJ 

z 


(/) 

Ui 

UJ 

c 

o  o 

UJ  o 

o  o 
o 
o  • 
o 

(O 

o  .. 

UJ 

z 


(/>  • 

CC  «  J 

o 

li.  UJ 

X  (/) 
</)  D  - 
^  K 

<o  <  X 

z  -so 

O  UI  u. 

^  I  a. 

K  S  -• 

<  UJ 
3  ♦ 

a 

UJ  X  J  u. 

<  o 

a  (N  u 

<  O  M  Ui 

Ui  •  >-  3 

Z  J 
^  X  < 
-I  «J  > 


UI 

z 


X 

o 


fU 

o 


UJ 

z 


o 

z 


H  u 


If 

H  > 


Ui  UJ 


tA 

^  to 

O  X 

z  •  o 

O  u. 

<  r 

3  ♦  U 

a 

UJ  X  U. 
O 

X  o 

<  o  UJ 
UJ  '  3 
Z  ^ 
—  X 
J  > 

UJ  II 

^ 


fU 

o 


(O 

UJ 

z 


UI 

z 


to 
UJ 
UJ 
X 

o  o 

UJ  O 

o  o 
o 
o  • 

o 

cn 

o)  .. 

M«  UJ 

^  z 

iO  M 

X  ^ 

o 

U.  UJ  K 
X  to 
to  3  — 
^  ^ 

-  <  X 
Z  •  X  o 
O  <*)  UI  u. 
M  I  a 
H*  S  "• 

<  UI  X 
3  ♦  H  U 

a 

UJ  X  3  U. 

<  o 

X  X  u 

<  O  UJ 

UJ  •  3 

Z  M  u 
*-•  X  < 
-I  u  > 

ui  If  UJ  UJ 

r  r  X 

I-  >  I-  »- 


Cl 

o 


UI 

z 

M  UI 

-I  z 

X  -J 

o 

u,  O 

z 

to  M 

o 

O  X 

z  •  o 

O  U. 

K  « 

<  X 

3  >  U 

a 

UJ  X  u. 

o 

X  o 

<  O  UI 
UJ  •  3 
Z  ^ 

< 

^  > 

UJ  II  UJ 

X  X 

H  >  I- 


APPENDIX  D 


This  appendix  contains  the 
temperature-  data  for  the  output 
The  data  was  hand-digitized  from 
published  source  14]. 


tosistance  versus 
contained  in  Appendix 
a  previously 
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temperature  (K)  resistance 
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